Spin polarization phenomena in dense neutron matter at a strong magnetic field 
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Spin polarized states in neutron matter at strong magnetic fields up to 10 18 G are con- 
sidered in the model with the Skyrme effective interaction. Analyzing the self-consistent 
equations at zero temperature, it is shown that a thcrmodynamically stable branch of solu- 
tions for the spin polarization parameter as a function of density corresponds to the negative 
spin polarization when the majority of neutron spins are oriented oppositely to the direction 
of the magnetic field. Besides, it is found that in a strong magnetic field the state with the 
positive spin polarization can be realized as a metastable state at the high density region in 
neutron matter. At finite temperature, the entropy of the thermodynamically stable branch 
demonstrates the unusual behavior being larger than that for the nonpolarized state (at 
vanishing magnetic field) above certain critical density which is caused by the dependence 
of the entropy on the effective masses of neutrons in a spin polarized state. 
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I. INTRODUCTION. BASIC EQUATIONS 

Neutron stars observed in nature are magne- 
tized objects with the magnetic field strength at 
the surface in the range 10 9 -10 13 G [l|. For a spe- 
cial class of neutron stars such as soft gamma-ray 
repeaters and anomalous X-ray pulsars, the field 
strength can be much larger and is estimated to 
be about 10 14 -10 15 G [2]. These strongly magne- 
tized objects are called magnetars 0] and com- 
prise about 10% of the whole population of neu- 
tron stars 0| . However, in the interior of a mag- 
netar the magnetic field strength may be even 
larger, reaching the values about 10 18 G [1, 
Under such circumstances, the issue of interest is 
the behavior of a neutron star matter in a strong 
magnetic field 041] ■ Further we will approximate 
the neutron star matter by pure neutron matter 
as was done, e.g., in the recent study [8|. As a 
framework for consideration, we choose a Fermi 
liquid app roach for description of nuclear mat- 
ter H, E3 and as a potential of nucleon-nucleon 
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interaction, we utilize the Skyrme effective forces. 

The normal (nonsuperfluid) states of neu- 
tron matter are described by the normal distri- 
bution function of neutrons f K1K2 = Tr ga^ 2 a K1 , 
where k = (p,cr), p is momentum, a is the pro- 
jection of spin on the third axis, and g is the den- 
sity matrix of the system 
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be assumed that the third axis is directed along 
the external magnetic field H. The self-consistent 
matrix equation for determining the distribution 
function / follows from the minimum condition 
of the thermodynamic potential 0] and is 



/ = {exp(Yoe + Y 4 ) + I}" 1 = {exp(Yo£) + I}" 1 •(!) 



Here the single particle energy e and the quantity 
I4 are matrices in the space of k variables, with 

*4kik2 = Y^kiks, Y o = l/ T , and *4 = -fJ'o/T be- 
ing the Lagrange multipliers, [j,q being the chem- 
ical potential of neutrons, and T the tempera- 
ture. Given the possibility for alignment of neu- 
tron spins along or oppositely to the magnetic 
field H, the normal distribution function of neu- 
trons and single particle energy can be expanded 
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in the Pauli matrices <7j in spin space 

/(p) = /o(pVo + /s(pV3, (2) 
e(p) = £o(p)o" + £3(p)o"3- 

Using Eqs. (pQ) and ([2D, one can express ev- 
idently the distribution functions /o, /3 in terms 
of the quantities e: 

/o = ^{n(uj + ) + n(w_)}, (3) 
/ 3 = ^{n(uj + ) - n{u)J)}. 

Here n(uj) = {exp(Yo^) + 1} and 

w± = fo±6, (4) 

Co = ^0 - MO) ^3 = £3- 

As follows from the structure of the distribu- 
tion functions /, the quantities uj± play the role 
of the quasiparticle spectrum and correspond to 
neutrons with spin up and spin down. The dis- 
tribution functions / should satisfy the norma- 
lization conditions 

££/o(p) = ft (5) 
p 

y Yl /s(p) = et~Q^ = A e- ( 6 ) 
p 

Here g = g^ + q± is the total density of neutron 
matter, g^ and g± are the neutron number den- 
sities with spin up and spin down, respectively. 
The quantity Ag may be regarded as the neutron 
spin order parameter. It determines the magne- 
tization of the system M = /x n Ag, fj, n being the 
neutron magnetic moment. The magnetization 
may contribute to the internal magnetic field B = 
H -\- AkM . However, we will assume, analogously 
to Refs. [a,l8|], that the contribution of the magne- 
tization to the magnetic field B remains small for 
all relevant densities and magnetic field strengths, 
and, hence, B ~ H. This assumption holds true 
due to the tiny value of the neutron magnetic mo- 
ment fin = -1.9130427(5)^ « -6.031 • 10~ 18 
MeV/G [14] (/xjv being the nuclear magneton) 
and is confirmed numerically in a subsequent in- 
tegration of the self-consistent equations. 

In order to get the self-consistent equations 
for the components of the single particle energy, 



one has to set the energy functional of the system. 
In view of the above approximation, it reads [12] 

E(f) = E (f, H) + E int (f) + E field , (7) 

E (f, H) = 2Y, eo(p)/o(p) - 2/%# £ ^(p)> 
p p 

Eint(f) = ^{£o(p)/o(p) + e 3 (p)/ 3 (p)}, 
p 

E - H \ 

Efield ~ g^V, 

where 

£o(p) = ^£fl?(k)/o(q), k= (8) 
q 

q 

2 

Here £o(p) = is the free single particle 

spectrum, mo is the bare mass of a neutron, 
Uq(\c), Ui(k) are the normal Fermi liquid (FL) 
amplitudes, and £o>£3 are the FL corrections to 
the free single particle spectrum. Note that in 
this study we will not be interested in the total 
energy density and pressure in the interior of a 
neutron star. By this reason, the field contribu- 
tion Ef ie ici, being the energy of the magnetic field 
in the absence of matter, can be omitted. Using 
Eq. ([7]), one can get the self-consistent equations 
in the form [13] 

£o(p) = £o(p) + £o(p) - Mo, (9) 
6(p) = -fJ-nH + e 3 (p). 

To obtain numerical results, we utilize the 
effective Skyrme interaction. The normal FL am- 
plitudes can be expressed in terms of the Skyrme 
force parameters @, Gi| : 

C^(k) = 2t (l-x ) + |^(l-x 3 ) (10) 

+ -^[t 1 (l-xi) + 3t 2 (l + x 2 )]k 2 , 
t^(k) = -2t (l - x ) - |/(1 - x a ) + p- 

(11) 

• [t 2 (l + x 2 ) - *i(l - xi)]k 2 = a n + 6 n k 2 . 
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Further we do not take into account the effective 
tensor forces, which lead to coupling of the mo- 
mentum and spin degrees of freedom, and, corre- 
spondingly, to anisotropy in the momentum de- 
pendence of FL amplitudes with respect to the 
spin quantization axis. Then 



(12) 



2m n 

^ = - fin H+(a n + b n ^-)=^ + ^(^) 3 , (13) 



p 2 , Ag i b n 2 



4 ' 4 16 
where the effective neutron mass m n reads 

h 2 h 2 ■ 



2m r 



2l7lQ 



+ f[ti(l-x 1 ) + 3t 2 (l + x 2 )], 



(14) 



and the renormalized chemical potential fj, should 
be determined from Eq. ([5]). The quantity (q 2 )3 
in Eq. (|13|) is the second order moment of the 
distribution function fy: 



V„2 

q 



(15) 



In view of Eqs. (fl~2]) . (fT3|) . the branches u± = uj a 
of the quasiparticle spectrum in Eq. ([!]) read 



2m a 



fi + a(-fi n H + 



a n Ag b n 2 , 
~ + 16 (q ^ 



(16) 

where m a is the effective mass of a neutron with 
spin up (a = +1) and spin down (er = —1) 



h 2 h 2 Pa Q fj 



(17) 



• [ti(i - xi) + t 2 (i + x 2 )], = eta)- 

Thus, with account of expressions ([3]) for 
the distribution functions /, we obtain the self- 
consistent equations ©, ([6]), and (fT5]) for the ef- 
fective chemical potential fi, spin order param- 
eter Ag, and second order moment (q 2 )3- To 
check the thermodynamic stability of different so- 
lutions of the self-consistent equations, it is nec- 
essary to compare the corresponding free energies 
F = E — TS, where the entropy reads 



S = - 2^ 2^ {n(0 In n(u a ) (18) 
p o-=t,4 

+ n{u a ) In n(w ff )}, n(uj) = 1 — n{u). 



II. ANALYSIS OF THE 
SELF-CONSISTENT EQUATIONS 

In solving numerically the self-consistent 
equations, we utilize SLy7 Skyrme force [l5|], con- 
strained originally to reproduce the results of mi- 
croscopic neutron matter calculations. We con- 
sider magnetic fields up to the values allowed by 
the scalar virial theorem. For a neutron star with 
the mass M and radius R, equating the mag- 
netic field energy E H ~ (AirR 3 /3)(H 2 /8vr) with 
the gravitational binding energy Eq ~ GM 2 /R, 
one gets the estimate H max ~ ^■(6G) 1 / 2 . For 
a typical neutron star with M = 1.5M Q and 
R = lO~ 5 i?0, this yields for the maximum value 
of the magnetic field strength H max ~ 10 18 G. 
This magnitude can be expected in the interior 
of a magnetar while recent observations report 
the surface values up to H ~ 10 15 G [l||. 

Fig. [2 shows the neutron spin polarization 
parameter II = Ag/g as a function of density for 
a set of fixed values of the magnetic field at zero 
temperature. At H = 0, the self-consistent equa- 
tions are invariant with respect to the global flip 
of neutron spins and we have two branches of so- 
lutions for the spin polarization parameter, n,|(£>) 
(upper) and ITq (g) (lower) which differ only by 
sign, 11+ (g) = -n (g). At H ^ 0, the self- 
consistent equations lose the invariance with re- 
spect to the global flip of the spins and, as a con- 
sequence, the branches of spontaneous polariza- 
tion are modified differently by the magnetic field. 
The lower branch Ui(g), corresponding to the 
negative spin polarization, extends down to the 
very low densities. There are three characteristic 
density domains for this branch. At low densities 
Q ^ 0.5^0; the absolute value of the spin polariza- 
tion parameter increases with decreasing density. 
At intermediate densities 0.5 go ^ Q ^ 3go> there 
is a plateau in the Ili(g) dependence, whose char- 
acteristic value depends on H, e.g., IIi ?» —0.08 
at H = 10 18 G. At densities g > 3go, the magni- 
tude of the spin polarization parameter increases 
with density, and neutrons become totally polar- 
ized at g 6 go- 

It is seen also from Fig. [T] that beginning 
from some threshold density the self-consistent 
equations at a given density have two positive so- 
lutions for the spin polarization parameter (apart 
from one negative solution). These solutions be- 
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FIG. 1: (Color online) Neutron spin polarization pa- 
rameter as a function of density at T = and dif- 
ferent magnetic field strengths for SLy7 interaction. 
The branches of spontaneous polarization IT^~ , II J are 
shown by solid curves. 



FIG. 2: (Color online) Energy per neutron as a func- 
tion of density at T = for different branches Ui(g)- 
113(g) °f solutions of the self-consistent equations at 
H = 10 18 G, including a spontaneously polarized 
state. 



long to two branches, 112(5) and 113(g), charac- 
terized by the different dependence from density. 
For the branch 112(5), the spin polarization pa- 
rameter decreases with density and tends to zero 
value while for the branch 113(5) it increases with 
density and is saturated. These branches appear 
step-wise at the same threshold density g t h de- 
pendent on the magnetic field and being larger 
than the critical density of spontaneous spin in- 
stability in neutron matter. For example, for 
SLy7 interaction, g t h ~ 3.80 go at H = 5 ■ 10 17 G, 
and ^th ~ 3.92 qq at H = 10 18 G. The magnetic 
field, due to the negative value of the neutron 
magnetic moment, tends to orient the neutron 
spins oppositely to the magnetic field direction. 
As a result, the spin polarization parameter for 
the branches 112(5), 1^3(5) with the positive spin 
polarization is smaller than that for the branch of 
spontaneous polarization FIq , and, vice versa, the 
magnitude of the spin polarization parameter for 
the branch IIi(5) with the negative spin polar- 
ization is larger than the corresponding value for 
the branch of spontaneous polarization I1q . Note 
that the impact of even such strong magnetic field 
as H = 10 17 G is small: The spin polarization pa- 
rameter for all three branches IIi(5)-ri3(5) is ei- 
ther close to zero, or close to its value in the state 
with spontaneous polarization, which is governed 



by the spin-dependent medium correlations. 

Thus, at densities larger than 5th, we have 
three branches of solutions: one of them, IIi(5), 
with the negative spin polarization and two oth- 
ers, 112(5) an d 113(5), with the positive polariza- 
tion. In order to clarify, which branch is thermo- 
dynamically preferable, one should compare the 
corresponding free energies. Fig. [2] shows the en- 
ergy per neutron function of density at T = 
and H = 10 18 G for these three branches, com- 
pared with the energy per neutron for a sponta- 
neously polarized state [the branches U^(g)}. It 
is seen that the state with the majority of neutron 
spins oriented oppositely to the direction of the 
magnetic field [the branch 111(5)] h as a lowest en- 
ergy. However, the state, described by the branch 
113(5) with the positive spin polarization, has the 
energy very close to that of the thermodynami- 
cally stable state. This means that despite the 
presence of a strong magnetic field H ~ 10 18 G, 
the state with the majority of neutron spins di- 
rected along the magnetic field can be realized as 
a metastable state in the dense core of a neu- 
tron star in the model consideration with the 
Skyrme effective interaction. Note here that in 
the study [1] of neutron matter at a strong mag- 
netic field only thermodynamically stable branch 
of solutions for the spin polarization parameter 
was found in the model with the SLy7 Skyrme 
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FIG. 3: (Color online) The entropy per neutron for 
the III branch, measured from its value in the nonpo- 
larized state (at H — 0), as a function of density at 
H = 10 18 G and different temperatures. 

interaction. 

One can consider also finite temperature ef- 
fects on spin polarized states in neutron matter at 
a strong magnetic field. Calculations show that 
the influence of finite temperatures on spin polar- 
ization remains moderate in the Skyrme model, 
at least, for temperatures relevant for protoneu- 
tron stars (up to 60 MeV). An unexpected mo- 
ment appears when we consider the behavior of 
the entropy of spin polarized state function 



of density. Fig.[3]shows the density dependence of 
the difference between the entropies per neutron 
of the polarized (the ITi branch) and nonpolarized 
(at H = 0) states at different fixed temperatures. 
It is seen that with increasing density the differ- 
ence of the entropies becomes positive. It looks 
like the polarized state in a strong magnetic field 
beginning from some critical density g s is less or- 
dered than the nonpolarized state. Such unusual 
behavior of the entropy was found also in the ear- 
lier works for spontaneo usly polarized states in 
neutron [l7] and nuclear (la . Qj| matter with the 
Skyrme and Gogny effective forces, respectively. 
Providing the low temperature expansion for the 
entropy in Eq. (|18p . one can get the condition 
for the difference between the entropies per neu- 
tron of the polarized and nonpolarized states to 
be negative in the form 

^t(i + n)l + ^(i-n)l-2<o. (19) 

For low temperatures, it can be checked nu- 
merically that this condition is violated for the ITi 
branch of the spin polarization parameter above 
the critical density g s being weakly dependent on 
temperature. 
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